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We study the effect of Faraday rotation due to a uniform magnetic field on the polarization of 
the cosmic microwave background (CMB). Scalar fluctuations give rise only to parity-even E-type 

■ polarization of the CMB. However in the presence of a magnetic field, a non-vanishing parity-odd 
' B-type polarization component is produced through Faraday rotation. We derive the exact solution 

f"**) , for the E and B modes generated by scalar perturbations including the Faraday rotation effect of 

■ a uniform magnetic field, and evaluate their cross-correlations with temperature anisotropies. We 
j 1 ' compute the angular autocorrelation function of the B-modes in the limit that the Faraday rotation 

I is small. We find that uniform primordial magnetic fields of present strength around Bo = 10 _9 G 

. rotate E-modes into B-modes with amplitude comparable to those due to the weak gravitational 

lensing effect at frequencies around v = 30 GHz. The strength of B-modes produced by Faraday 
rotation scales as Bq/v 2 . We evaluate also the depolarizing effect of Faraday rotation upon the cross 
correlation between temperature anisotropy and E-type polarization. 

PACS numbers: 98.70.Vc,98.80.Es,98.80.Cq 
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I. INTRODUCTION 
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^ ■ 

The detection of a parity-odd component of the CMB polarization, called B-mode, can provide a unique tool to 
probe sub-dominant sources of CMB perturbations, which are primarily due to scalar energy-density fluctuations. 
Scalar perturbations do not produce B polarization, which is only excited by either tensor or vector modes Q, . It 
is thus important to identify all potential sources of B polarization. A significant background of B-modes comes from 
the effect of gravitational lensing on the CMB by the matter distribution, which transforms E into B polarization 
Cosmological gravitational waves, as those predicted by inflationary models of the early universe, would also 
imprint B-polarization upon the CMB. It has been pointed out that the signature of an inflationary gravitational- 
wave background can only be detected by B-polarization measurements if the tensor to scalar ratio r > 10~ 4 , which 
corresponds to an energy scale of inflation larger than 3 x 10 15 GeV 0,01 , because otherwise their effect would be 
obscured by the weak lensing background. Quite recently, however, a better technique to clean polarization maps from 
the lensingeffect has been proposed, which would allow tensor-to-scalar ratios as low as 10 -6 , or even smaller, to be 
probed |a,l3- Other secondary contributions to the B-type polarization arising during the reionization stage, though 
of much smaller amplitude, have been considered in [g. B-modes of polarization are also produced by secondary 
vector and tensor modes generated by the non- linear evolution of primordial scalar perturbations 0- Primordial 
magnetic fields can also generate B-mode polarization upon the CMB, due to the vorticity they induce upon the 
photon-baryon fluid O 111 [13 • 

We will consider here another source of B-type polarization, that is the Faraday rotation in a magnetic field of 
initial E-type polarization. If a primordial magnetic field is present at the time of the last scattering of the microwave 
photons, it causes Faraday rotation of the direction of the linear polarization which generates B-modes from 
initial E-modes of polarization |l4j. In Ref. 0] a cross-correlation between the parity-odd B-modes generated by 
Faraday rotation and the parity-even temperature anisotropies O was estimated under simplifying approximations. 
0-B or E-B correlations also exist if there are other parity violating processes, for instance in the case in which the 
vorticity in the photon-baryon fluid is induced by helical magnetic fields generated by parity-violating electroweak 
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interactions |15l llq , or if parity- violating interactions additional to those in the standard electroweak theory exist 

E3 

Large scale magnetic fields have been observed in galaxies, clusters and superclusters, but their origin is still unknown 
0GJ1. They may have formed from dynamo amplification of a small initial seed field poj , or alternatively they can 
arise from the adiabatic compression during structure formation of a cosmological field of order 10~ 10 — 10~ 9 G today 
|2*lj . The existence of these very large scale cosmological fields, although is not excluded by present observations, is 
difficult to be determined observationally. 

Here we derive the exact solution for the E and B modes generated by scalar perturbations including the Faraday 
rotation effect of a uniform primordial magnetic field, and discuss their properties. We evaluate the effect of Faraday 
rotation upon the cross correlation O-E and the induced correlation 0-B. We then calculate, to first order in the 
Faraday rotation effect, the angular correlation function of the B-modes induced by Faraday rotation of scalar E- 
modes of polarization in a standard cosmological model, and compare it with that due to other potential sources of 
B-modes, namely weak gravitational lensing, gravitational waves, and the effect of Faraday rotation due to magnetic 
fields in clusters [HEl- 



II. FARADAY ROTATION OF CMB POLARIZATION 



Polarization of the CMB arises from Thomson scattering of anisotropic radiation by free electrons. One way of 
describing linear polarization is by means of the Stokes parameters Q and U [24j . A magnetic field induces a rotation 
of the polarization plane of a linearly polarized wave as it propagates through an ionized medium. This effect is known 
as Faraday rotation. If there were a primordial magnetic field at the recombination epoch, or after the Universe was 
reionized, it would induce Faraday rotation on the CMB photons, mixing U and Q modes of polarization. The 
evolution equations for their Fourier modes read [T3L I25I |26| 

3 

Q + ikfiQ = -tQ + -f(l - t i 2 )S p + 2uj B U 

U + ikfxU = -tU-2lu b Q, (1) 

where n = n ■ k is the cosine of the angle between the photon direction and the Fourier wave-vector, and f is the 
differential optical depth (r = naxa/ao with <jt the Thomson scattering cross section, n the free electron density, 
and the derivative is taken with respect to conformal time dr\ — dt/a(t)). We define lob = 3£? • nc 2 f/167r 2 i/ 2 e as the 
Faraday rotation rate of radiation with frequency v in a magnetic field B (c is the speed of light and e the electron 
charge). Polarization is sourced by quadrupole anisotropies in the temperature fluctuations 0. S p in eq. Q is the 
source term, which reads S p = (02 + C?2)/5 + Qo if and Q are expanded in terms of Legendre polinomials as 
= J2e(~ iYQiPtif*) ( ancl a similar expansion for Q)[3(J. 

The equations above have the following formal integral solutions 

q rvo 

{Q±iU){k, m ) = -(l-iJ?) dr)e- ik ^°-^Te- T S p (ri)e* iFTB -*. (2) 
4 Jo 

Here t(jj) = J vo drj'f(r]') is the optical depth between time 77 and present time 770, and the parameter F is defined 

as 



3 Be 2 „ / B \ /10GHz x 2 



F= 8^~°- 7 U=^J {—) ' (3) 

F gives a measure of the typical Faraday rotation between collisions. Notice that the strength of a primordial 
magnetic field is expected to evolve in time cx a~ 2 due to flux conservation. Since v cx a -1 , then F (expressed here 
numerically in terms of the present values of B and v) does not change with time. 



III. E AND B MODES 



The generation of temperature andpolarization anisotropies in the CMB from gravitational perturbations has been 
studied in detail by several authors 0, H, 123, |28| . A simple and powerful formalism is the total angular momentum 
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method [271 l28| . which includes the effect of scalar, vector and tensor modes on an equal footing. We are going to 
extensively use the results and notation of Ref. ^2^ . 

The temperature and polarization fluctuations are expanded in normal modes that take into account the dependence 
on both the angular direction of photon propagation n and the spatial position x as 



(27T) 

(Q±iU)(r,,x,n) = [ ± <^ m) ) ia^T > 



( 4 ) 



with spin s = describing the temperature fluctuation and s = ±2 describing the polarization tensor and to 
0, ±f,±2 corresponding to scalar, vector and tensor perturbations, respectively. E^ and are the angular 

moments of the electric and magnetic polarization components and 



s GT(x,h) = Hr^^ m (»K k , (5) 

is the basis of the expansion, in terms of the spin-weighted spherical harmonics. 

The set of Boltzmann equations for the evolution of , E^ 1 ' and B^ , as well as their integral solutions are 
deduced in Ref. [27j. In the presence of a uniform magnetic field, Faraday rotation modifies the Boltzmann equations 

satisfied by E^ and B^ m \ We quote in the Appendix the modified equations for completeness. We can nevertheless 
sidestep them inverting equation (@J to derive an integral solution for E^™ 1 ' ± iB^ from the integral solution for 
Q ± iU in eq. 10) 



where the source term S p has been expressed in terms of p(°) = (6^ 0) - V&E { ° ] )/\Q = S p /2 [27|. First the exponential 
with the magnetic field is expanded in terms of spherical Bessel functions and ordinary spherical harmonics as 

eTl FrB-n = ^^(f^y™ (n)Y ( m * (B) . (7) 



Lm 



Then the product of one ordinary and one spin- weighted spherical harmonic is written as a sum of terms proportional 
to just one s YI n {n) using the Clebsh-Gordan relations for the addition of angular momentum |2*7| : 



• vmiW V ™,2\ V (2^1 + 1)(2£ 2 + 1) ST^ 16 \D I) \ 

,sx Y ii ) \s2 Y i 2 ) = 4^ (h,h;rni,m 2 \h,t2;i,Tn) 



The term exp (— ifc/i(ryo — T])) is also expanded in terms of spherical Bessel functions, with the £3 axis chosen in the 
direction of k (which leaves only Y^(n) in the expansion) and the angular integration is performed. The final result is 

E^ m) ±iB} m) = (-l) m+1 V24^(2£+ 1) / d V Te- T P {Q) J2( i y^'' Te ' Y e r (B)ji{FT)VWTTx 

Jo t , t „ 



x <£, £'; -to, m\£, if; £",0><£, £'; ±2, 0\£, £'; £", ±2 > e%> (k( m - 77)) . (8) 

where ef\x) = y/3(l + 2)!/8(Z - 2)\j e (x)/x 2 . 

Using the relation < £, £'; -2,0\£, £'; £" , -2 >= (-if-W < £,£';+2,Q\£, £';£", +2 > between Clebsh-Gordan 
coefficients, individual expressions for E^ and iB/j immediately follow. They are similar to the right hand side of 
eq. © choosing the upper sign, with an extra factor (1 + (— l) e ~ e )/2 in the case of E^ 1 ' and (1 — ~ e )/2 in 

the case of iBi m ^ inside the sum. 
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Notice that in the absence of a magnetic field the standard solution for E^ is recovered, and also B^ = 

as it should. Instead, in the limit of relatively large magnetic field (F >> 1), both and B^ tend to zero due 

to the rapidly oscillatory nature of the spherical Bessel functions. This is a manifestation of the depolarizing effect of 
differential Faraday rotation [2{| . 

The presence of the magnetic field, which breaks the invariance under rotations, gives rise to a coupling between 
different m multipoles. This is evident from eq. JSJ, where a full range of m polarization multipolcs is obtained from 
the initial scalar m — modes, as well as from the Boltzmann equations displayed in the Appendix. 



IV. FARADAY ROTATION EFFECT UPON CROSS-CORRELATIONS 

In order to calculate correlation functions one has to integrate over wavevectors k, for which our expressions for 
e[ m \ E^ m) and B { ™\ valid in a frame with £3 || fc, need to be rotated to a frame with fixed axis. In the absence of 
a magnetic field, statistical homogeneity and isotropy makes the angular integrations simple, and the angular power 
spectrum of the correlation functions can be calculated as 

with X, Z — O, E, B. Statistical homogeneity and isotropy is reflected in the fact that 

< Xi m Z^ rn i >= S mm 'Su'Cf z , (10) 
where Xi m are the coefficients of the multipole expansions in the fixed basis: 

(Q±iU)(n) = Ee,m( E em ± %B tm ) ±2 Y e m . 

A uniform magnetic field breaks invariance under rotations and defines a preferred direction in space. We can take 
the direction of the magnetic field as the S3 direction of the fixed basis. The multipole coefficients Ei m ± iBi rn in this 
basis can be calculated through a rotation of those given by eq. © , as 

E lm ± iB lm = J 0e*»£(f£ m ' ) ± iB^) -m'Y e m *(k) . (12) 

m' 

Here E { ™' ] ± iBf 1 ^ are given by eq. ((HJ) with B replaced by k in the argument of YIP , since now the angle between 
B and k is measured from the fixed axis along B. The factor summed over m! in eq. i|12|) corresponds to the 
rotation of the polarization multipoles corresponding to each wavevector k mode to the fixed axis system, multiplied 
by yj (2£ + l)/47r. The other factors arise from the definition of the basis ± 2 GJ l used in eq. (QJ). The product of 
two spherical harmonics appearing after replacing eq. JSJ in this expression can be expanded with the relations for 
addition of angular momentum, with the help of the relation (— l) m m Y®(8, <fi) = Yf 1 ^, 0), and the sum over m! can 
be done using the orthogonality properties of the Clebsch-Gordan coefficients. The result can finally be expressed in 
terms of the analogous multipole coefficients in the absence of a magnetic field as [37j 



Elm ± iBh n — 



■ir-' + l ± (-l)' 



R{£m,Z',£"m)E { p m . (13) 



Here Eg m coincides with the expression for the multipole coefficients due to scalar fluctuations in the absence of a 
magnetic field, except for an extra term (2f + l)j>(i<Y) inside the time-integral that defines it. Explicitly: 

C ) = -^^TT / P^e rkS Yr(k)V6(2£+l) f° drp-e^W + l)M^)P<°>fa)e<°>(fcfa, - r?)) . (14) 







lm " 2£+l J 2tt 3 
The coefficients R(£m, £', £"m") are given, in terms of Clebsch-Gordan coefficients, by 



R{£m, £', £"m") = / ^ j== < t'\ m - m " \ m"\£'£"; £, m >< £', £"; 0, 2\£'£"; £,2 > . (15) 
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FIG. 1: Temperature-polarization (E-mode) cross-correlation in a ACDM background cosmology in the case of no magnetic 
field (thin solid line) and in the case of Faraday rotation in a uniform magnetic field with F = 1 (thick solid line) and F = 3 
(dashed line). F is defined in eq. ®- For large F the depolarizing effect of differential Faraday rotation is quite noticeable. 



Correlation functions including the Faraday rotation effect can now be immediately calculated in terms of their 
values in the absence of the magnetic field, as given by eq. QUI II . The multipoles of the cross-correlation between 
temperature anisotropies and E and B modes can be written as 



(-If 



- 1 



(16) 



This expression is strictly valid in a frame with 63 || B, and the fact that only objects with the same m are correlated 
reflects the cylindrical symmetry of the problem around this axis. The cross-correlation in any other system can 
be obtained by noticing that the multipoles in a new system are related to the ones in the £3 || B system by 
X' em — ^jAi:/2£ + 1 J2m' Xtm' - m 'Yi m *{B). In the general system the cross correlation results 



)> = 



-iy 



1 ± (-1)' 



1 



4tt 



2C 



1 



Y " 
1 fi 



\B)Cl 



9E l 



Thus, in a general system the correlations between different m values do not vanish. This fact can be used to identify 
the presence of a symmetry axis, that in this case corresponds to the direction of the magnetic field. 

It is clear from the expressions in eqs. ()16|) and (|17J) that Faraday rotation in a uniform magnetic field generates 
non-diagonal (£\ 7^ £2) correlations 0-B out from initially diagonal correlations 9-E. The difference between £\ and 
£2 must be odd, and contributions to the case \i 2 — £i\ — 2k + 1 are increasingly suppressed for larger k. They are in 
fact of order larger or equal than 2fe + 1 in the parameter F. This is because the coefficient R vanishes if |£a — £\\ > £' 
and since (F) oc F e for small F. Likewise, Faraday rotation induces E-B correlations with £ 2 — £\ odd. It also 
distorts the diagonal 0-E and E-E correlations, and induces non-diagonal terms with \l% — 1\\ = 2k, proportional at 
least to F 2k for small F. 

In the system with £3 || B, we can actually calculate to all orders the effect of Faraday rotation upon the diagonal 
part of the correlation 0-E. Indeed, since J2 m R{^ m i 1 = (2^ + l)^'o we find that 



c: 



0E — 



1 

2£TT 



2J < &£ m E em >— c, 



BE 11 



(18) 



In other words, the diagonal terms in the cross-correlation 0-E can be calculated, to all orders, exactly as in the 
case of no magnetic field with the addition of a term jo(FT) in the integrand along the line of sight for Eg m . Notice 
however that this equivalence is not valid for individual contributions to Ci with different values of m, since they also 
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FIG. 2: Visibility function g(rf) (solid line) as a function of conformal time in the background ACDM cosmology, calculated 
with CMBFAST, and multiplied by jo(Fr) with F — 3 (dotted line), and multiplied by r (dashed line). The left panel displays 
times around hydrogen recombination and photon decoupling. The right panel displays times after reionization, which was 
assumed to take place at an optical depth r re ion = 0.17. 



get additional contributions from higher E) 1 which, however, add to zero in the sum over m. We can readily evaluate 
expression (|18l) numerically using CMBFAST 29), account taken of the extra term jo(Fr). As background cosmology 
a spatially-flat A-CDM model was assumed, with 70% of the energy density in the form of a cosmological constant, 
25.6% in cold dark matter, 4.4% in baryons, a Hubble constant Ho = 71 km/sec Mpc, adiabatic scale invariant 
scalar fluctuations, and an optical depth since reionization r re i 0n = 0.17, as suggested by WMAP measurements of 
temperature-polarization correlations [3lJ . The result is plotted in Figure ^ m the standard case of no Faraday 
rotation (thin solid line) , and in the case of Faraday rotation factors F = 1 (thick solid line) and F = 3 (dashed line) . 

Notice that there is no linear effect of Faraday rotation upon Cf E for small F, the lowest order correction being 
quadratic and with negative sign. For large F the oscillatory nature of jo (-FY) leads to a significant reduction in 
the amplitude of the cross-correlation, reflecting the depolarizing effect of large differential Faraday rotation [25| . A 
shift towards larger angular scales (smaller I) in the positions of the peaks is also noticeable in Fig. ^ This can be 

interpreted as the consequence of the fact that the integrand in eq. (|14f> for E^ contains the factor g(j])jo(Fr), where 
17(77) = rexp (— t) is the visibility function (the probability that a photon last-scattered within drj of rf). The visibility 
function around the time of hydrogen recombination can be well approximated by a Gaussian centered at the time 
of photon decoupling. The polarization produced by Thomson scattering at recombination can be well approximated 
as proportional to the time derivative of temperature inhomogeneities around decoupling and to the width of the 
visibility function. The additional factor jo(Fr) in the integrand of eq. (|14|) not only reduces the value of the integral, 
but also shifts to later times the location of the largest contribution. We display this effect in Fig. |2 where we plot 
the visibility function g(r) (solid line) and the factor g(r])jo{FT) for F — 3 (dotted line), both around the time of 
photon decoupling (left panel) and after reionization (right panel) . The shift to later times in the location of the peak 
in g(r])jo(FT) in the left panel translates into a shift towards larger angular scales in the position of the peaks in the 
correlation function. 

Faraday rotation has a relatively smaller effect upon E modes produced after reionization because r is not very 
large (smaller than 0.17 in our numerical examples), and thus jo(Fr) is closer to unity unless F is larger. 

The strong damping of 0E correlations that occurs when F is of order unity or larger also affects 0B, EB, EE 
and BB correlations. It is a consequence that differential Faraday rotation across the last scattering surface (or after 
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FIG. 3: Temperature-polarization (B-mode) cross-correlation (thick solid line) in a ACDM background cosmology due to 
Faraday rotation of E-modes, scaled by the factor F defined in eq. 10, assumed to be smaller than unity. See the text for the 
definition of C® B , which is actually a non-diagonal correlation. The correlation O-E in the absence of magnetic field is also 
displayed for reference (thin solid line). 



rcionization if F were large enough) has a net depolarizing effect. 

Corrections to Cf E smaller than those considered here exist due to the fact that Faraday rotation affects the 
quadrupole term in the temperature anisotropy source [2£j], an effect that we have neglected. 

It is clear from Fig. ^ that Faraday rotation factors of order unity would imprint a measurable signature upon the 
temperature-polarization cross-correlation, which is moreover frequency-dependent. Primordial magnetic fields with 
strength 10~ 8 G (scaled to present values) can easily be ruled out with measurements at frequencies of order 30 GHz. 



CORRELATION FUNCTIONS FOR SMALL FARADAY ROTATION 



In the (likely to be realistic) case of small Faraday rotation (F « 1), it suffices to take t' = 1 in cq. (|16(l or Ijl7|l . 
To this order, the cross-correlation 9-B in the system with <=3 || B reads 



9£' 



. v^i + i^-^V^i + i) 2 - 4 , - 

0£ 2 /i + l ~ : — ■- - r Of 



(19) 



(h + i)V2£i + iVW+3 ' ~"' 2 '" L £W2h ~ V2ii + 1 

In a general frame there are also correlations among multipoles with values of m differing by ±1. Notice that a similar 
expression holds for the cross correlation E-B, with C^ <0) ^ ( ' in place of Cf^ . 

In Fig. |3 we plot, as a measure of the strength of the cross-correlations 0-B, the quantity 



a 



6B _ ^ei (11 



2( 



+ i ^ 



2 — m 2 \[P 



-—I 



Isjli - ly/2t- 



(20) 



which represents an average over m of the amplitude of the non-diagonal correlations. We plot it scaled down by the 
factor F, assumed to be small. Notice that Cf coincides with the power spectrum for 0-E correlations produced 
by scalar fluctuations in the absence of a magnetic field except for an extra factor 3ji (Ft) ~ Ft in the time integral 
that leads to E. This can be interpreted as the fact that the B-modes are the result of the E-modes produced by 
Thomson scattering up to a time r), further Faraday rotated by a factor Ft since, and integrated over 77. Fig. [^displays 
the factors g(r]) and Tg(rj) in the background cosmology chosen for illustration, both around photon decoupling (left 
panel) and after reionization (rigth panel). rg(rj) has its peak shifted to earlier times around decoupling compared 
with the visibility function. This effect shifts the peaks in the correlation functions to smaller angular scales (larger 
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FIG. 4: Angular autocorrelation multipoles of B-modes of CMB polarization generated by Faraday rotation of E-modes 
generated by scalar perturbations in a ACDM cosmology, scaling away the factor F 2 (thick solid line). The autocorrelation of 
E-modes from which they originate is also displayed (thin solid line) for comparison purposes. Also shown are the autocorrelation 
of B-modes due to weak gravitational lensing (dashed line) and to gravitational waves in a model with tensor to scalar ratio 
r = 1CT 2 (dotted line). 

values of £). Faraday rotation effects after reionization are not only suppressed by a factor F, but also by the relatively 
smaller values of t, as displayed in the right panel. 

The diagonal contribution to the angular autocorrelation of B-modes is given, to lowest order in F and averaged 
over m, by 

„BB_ 1 „ D . o .1 (l+l) 2 -4 iMjdl 1 I 2 -4 gWgW 

Ll = 27TT < tm lm 3 {2i + 3)(£ + 1) Ce+1 + 3(2e + i)e Ce - 1 ■ {21 > 

We calculate it numerically using CMB FAST [29j, accounting for the extra factors Ft in the integrals along the 
line of sight for The result is displayed in Figure 0] for the same background A-CDM cosmological model as 

in the previous section. The quantity plotted (thick solid line) is 1(1 + \)Cf B /F 2 , with F given by eq. ©, so that 
the result should be scaled down by the appropriate value of F 2 , which we have assumed to be small. We also 
plot, for comparison purposes, the B-mode correlation due to weak gravitational lensing (dashed line), and that due 
to gravitational waves in a model with tensor to scalar ratio r — 1CP 2 (following the convention in e.g. [^2 for 
its definition). Also plotted in Fig. 0] to help estimate the strength of the effect, is 1(1 + l)Cf E due to the scalar 
fluctuations in the absence of Faraday rotation (thin solid line). 

The angular correlation multipoles of B-modes are roughly a factor F 2 /3 smaller than those of the E-modes from 
which they originate by Faraday rotation, at relatively large I. Notice also the shift in the peaks positions towards 
larger I. At small I the effect is relatively suppressed because the B-modes produced by Faraday rotation are smaller 
than the E-modes generated after reionization by an extra factor somewhat smaller than r r 2 eion (the optical depth to 
reionization) . 

It is clear from Fig.0]that Faraday rotation factors of order F = 0.1, which at v = 30 GHz correspond to a primordial 
magnetic field of present strength Bq = 10~ 9 G, induce B-modes with amplitude comparable to those generated by 
weak gravitational lensing. B-modes due to Faraday rotation in the magnetic fields of clusters of galaxies could have 
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a comparable strength (but different angular spectrum) around £ — 1000 only if i? c iustcrs ~ 2 /xG, uniform across the 
clusters pll^. 



VI. CONCLUSIONS 

The detection of polarization in the CMB [3(j and its correlation to the temperature anisotropies jjuj has become 
a very promising tool to test the early universe. Much observational effort is taking place to measure 0, E and B 
correlation functions. The cross-correlation function OE already gives information on the reionization history of the 
universe. The B-mode polarization is being looked upon specially as a way to probe a primordial gravitational wave 
background which would represent a clear signature of a period of inflation in the early Universe. Precision data 
on the CMB polarization is expected from the Planck Mission and future dedicated missions, such as NASA's 
Beyond Einstein Inflation Probe [34| or ground-based experiments, like BICEP [3^ and PolarBeaR. 

There are other processes beside a gravitational wave background that can originate B-mode polarization, more 
prominently the gravitational lensing conversion of a fraction of the dominant E-type into B-type polarization, which 
can however be significantly cleaned from the data allowing to probe inflationary models with tensor-to-scalar ratios 
r<lO- 6 0. 

We have studied here in detail the effect of the Faraday rotation produced by a large scale uniform primordial 
magnetic field upon the polarization of the CMB for the case of scalar density perturbations. The exact solution for 
the E and B modes of polarization in the presence of a uniform magnetic field is given in eq. (|13[) . Faraday rotation of 
the direction of the polarization generates a B-polarization component from the initial E-modes, which correlates with 
temperature anisotropies and E-polarization .14]. 0-B correlations are non-diagonal, and their strength is of order 
F times the original 0-E correlations, with the parameter F defined by eq. (J3J, as long as F is smaller than unity. 
The amplitude of the B-polarization angular autocorrelation function generated by Faraday rotation is comparable 
to the (uncleaned) weak lensing signal if F ~ 0.1, which corresponds to a magnetic field B ~ 10~ 9 G at a frequency 
v ~ 30GHz. Let us notice however, that the lensing signal, as well as that generated by gravitational waves, are 
independent of the frequency, while F scales as v~ 2 , thus multi- frequency measurements should easily distinguish the 
Faraday rotation contribution. 

Small Faraday rotation factors have a relatively small effect upon 0E and EE correlations, since they induce 
distortions roughly F 2 smaller than their values in the absence of a magnetic field. If F were of order unity or 
larger, all correlations 0E, 0B, EE, EB and BB are significantly damped, due to the additional depolarizing effect of 
differential Faraday rotation across the last scattering surface. The good agreement between WMAP measurements 
[31 1 of 0E correlations and its predicted value in a standard cosmology rules out values of F of order unity, which 
implies primordial magnetic fields of present strength Bq 10 _8 G at frequencies v =30 GHz. Future precision 
measurements should allow to probe smaller values of primordial magnetic fields through Faraday rotation effects. 



APPENDIX: BOLTZMANN EQUATIONS FOR E AND B 

We write down here the Boltzmann equations for the E and B-modes of polarization including the Faraday rotation 
effect of a constant magnetic field. We adopt the formalism of Ref. [2]} for the case without a magnetic field. The 
additional terms, that mix E and B modes, can be obtained by derivation of the defining equation Q, using eq. QJ, 
and the Clebsh-Gordan relations for products of spherical harmonics. It turns out that 
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where Ob and 0s are the polar and azimuthal angles respectively of the direction of the magnetic field in a system 
with e3 || k, and we have defined 
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